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Abstract: Longitudinal double-spin asymmetries of charged hadrons with high transverse 
momentum px have been measured in electroproduction using the Hermes detector at 
Hera. Processes involving gluons in the nucleon have been enhanced relative to others 
by selecting hadrons with pr typically above 1 GeV. In this kinematic domain the gluon 
polarization has been extracted in leading order making use of the model embedded in the 
Monte Carlo Generator Pythia 6.2. The gluon polarization obtained from single inclusive 
hadrons in the px range 1 GeV < px < 2.5 GeV using a deuterium target is (/i 2 )) = 

0.049 ± 0.034(star:) ± 0.010(sys-exp)tomt(sys-models) at a scale (^ 2 ) = 1.35 GeV 2 and 
(x) = 0.22. For different final states and kinematic domains, consistent values of — have 
been found within statistical uncertainties using hydrogen and deuterium targets. 
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1. Introduction 



In recent years a major goal in the study of Quantum Chromo-Dynamics (QCD) has been 
the detailed investigation of the spin structure of the nucleon and the determination of the 
partonic composition of its spin projection 

± = S S = ±.£Z + AG + Ll + L<;. (1.1) 

Here AS is the contribution of the quark and anti-quark helicities, AG is the contribu- 
tion of the gluon helicity, and L\ and are the quark and gluon orbital angular mo- 
menta, respectively, in a reference system where the nucleon has very large longitudinal 
momentum. The individual terms in the sum depend on the scale \j? and the renormal- 
ization scheme. Recent results from experiments |3| and fits in next-to-leading order 
(NLO) QCD ||, [5], [| |7| to helicity-dependent inclusive Deep-Inelastic Scattering (DIS) 



data [|, |, |, |, [HJ, y, [lj, [lj, [lj, £|, |T|, gg, gj| yield a value of AS ~ 0.2 - 0.4 at 
/x 2 = 4 GeV 2 in the MS scheme [19]. In contrast to the quark helicity distributions, the 
knowledge of the gluon helicity distribution function is still limited. There are no di- 
rect experimental determinations of parton orbital angular momenta. Most of the existing 
knowledge about AG(fj?) originates from next-to-leading order perturbative QCD (pQCD) 
fits to the helicity-dependent structure function g\{xs, Q 2 ) of the nucleon, where xb is the 
Bjorken scaling variable, which is in leading-order (LO) identified with the longitudinal 
parton momentum fraction x in the nucleon. In DIS the renormalization and factoriza- 
tion scales fi 2 are set equal to the photon virtuality Q 2 . Because the virtual photon does 
not couple directly to gluons (see Fig. [j]b), gi(xB,Q 2 ) is only weakly sensitive to gluons 
through the DGLAP evolution [p0| pij, |22| of the helicity-dependent Parton Distribution 
Functions (PDFs). At next-to- leading order pQCD, additional sensitivity to gluons arises 
from the Photon-Gluon Fusion (PGF) subprocess (see Fig. |l|b). However, the limitations 
on the precision and kinematic range in xg and Q 2 of the g\ measurements result in large 
experimental and theoretical uncertainties on the determination of the gluon helicity dis- 
tribution function Ag(x,fi 2 ). Results for AG(fi 2 ) = Jq Ag(x, fi 2 )dx from recent pQCD fits 
to inclusive DIS data M, ||, ||, |7j are typically of order 0.5 with uncertainties up to ±1. 

An alternative constraint on the extraction of AG(fj?) in NLO pQCD fits comes from 
the measurements of double-spin asymmetries in production of inclusive 7r° mesons or 
jets with high transverse momentum in polarized proton-(anti-)proton scattering. First 



measurements were performed by E704 |23[ and more recent data were obtained by Phenix 



\2§ and Star p5j at Rmc. The inclusion of the RHIC-data in recent NLO pQCD fits @ 
improves the accuracy on Ag significantly. One finds \Ag(x,Q 2 )\ smaller than 0.1, with a 
possible node in the distribution. This is driven mainly by the RHIC data, which constrain 
the magnitude of Ag(x) for 0.05 < x < 0.2, but cannot determine its sign as they mainly 
probe the product of the gluon helicity distribution at two x values. 

In order to increase the sensitivity to Ag(x,fi 2 ) in lepton- nucleon scattering, other 
observables besides the inclusive helicity-dependent structure function have been studied. 
These observables are expected to include a direct contribution from gluons. For example, 



-3- 



7 7 Q 7 £ 









Figure 1: Feynman diagrams for hard subprocesses: a) O{a s ) DIS, b) 0(oq) Photon-Gluon Fusion, 
and c) 0(a\) QCD Compton scattering. 



in hadron leptoproduction this gluonic contribution can be relatively enhanced by detecting 
charmed hadrons, or inclusive hadrons or hadron pairs at high transverse momenta px- 

Charmed hadron electroproduction is a suitable channel because it is dominated by the 
PGF subprocess [^6| and a hard scale is introduced by the mass of the charm-quark pair, 
which makes pQCD calculations of this process possible. For light final state quarks, the 
selection of hadrons with high px enhances the relative contribution of the gluon subpro- 



cesses and the relevant transverse momenta provide the scale (see Sect. 5.3). In the high 
Pt domain other calculable hard pQCD subprocesses such as QCD Compton (QCDC) 
scattering (see Fig. ||c) are relatively enhanced as well, whereas soft, non-perturbative pro- 
cesses are suppressed. Charm electroproduction is being investigated by Compass p6|, p7j . 
Inclusive single hadron leptoproduction was studied by E155 [pq| . Hadron-pair leptopro- 
duction at high p T was studied by Hermes J29|, SMC Q and Compass [3^]. For 
these experiments high px is in the range from one to a few GeV. 

Throughout this paper, the term "LO" is applied to all leading order subprocesses con- 
tributing to hadron production at nonzero px- These are the tree level processes at 0(a]), 
but also the quark scattering process j*q — > q (DIS) at 0(cP s ). While the former processes 
involve hard gluons, and can therefore involve substantial parton transverse momentum 
Pt in the hard scattering, in the latter process px is equal to zero, but hadrons acquire px 
from soft initial and final state radiation. This paper presents the LO extraction of the 
gluon polarization ^(%) from longitudinal double-spin asymmetries of charged inclusive 
hadrons measured in electroproduction using a deuterium target by Hermes at Hera. The 
contributions of signal and background have been determined by a Pythia Monte Carlo 
simulation, which includes LO pQCD as well as non-perturbative subprocesses. Consis- 
tency checks have been performed for different kinematic regions, different final states and 
using data from a hydrogen target. The data taken with the deuterium target correspond 
to an integrated luminosity three times larger than that taken with the hydrogen target, 
see table ffl. Compared to the previous Hermes publication [pjJ, which used measurements 
of hadron pairs of opposite charge from a hydrogen target, this analysis includes a much 
larger sample of single hadrons, and a significantly more comprehensive treatment of the 
underlying physics processes 33, 34]. 



This paper is organized as follows: In Sect. the experimental method is described, 
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in Sect. || the asymmetry results are given, in Sect. ^ the determination of ^ with a 
description of the physics model of the reactions is discussed, in Sect. || the Pythia Monte 
Carlo simulation is described, in Sect. ||] the determination of ^-(x,fi 2 ) is explained, and 
in Sect. ^ the summary and conclusions are given. 



2. The Hermes Experiment 



Positrons of momentum 27.6 GeV were stored in the Hera lepton ring at Desy. The 
initially unpolarized beam was transversely polarised by an asymmetry in the emission of 
synchrotron radiation associated with a spin flip (Sokolov-Ternov mechanism |35| ). The po- 
larization was rotated to the longitudinal direction for passage through an gaseous internal 
fixed target of longitudinally nuclear-polarized atoms. The scattered positron and hadrons 
produced were detected in a forward magnetic spectrometer. The beam helicity was re- 
versed periodically. The beam polarization was measured continuously by two independent 
polarimeters using Compton backscattering of circularly polarized laser light ]36|, |37[ . The 
average beam polarization for the data used in this analysis is shown in Tab. |l|. The tar- 



get [35] consisted of longitudinally nuclear-polarized pure atomic hydrogen or deuterium gas 
in an open-ended 40 cm long storage cell. The cell was fed by an atomic-beam source based 
on Stern-Gerlach separation combined with radio-frequency transitions of atomic hyperfine 
states |39| . The sign of the nuclear polarization of the atoms was chosen randomly every 
60 s (90 s) for the hydrogen (deuterium) target. The polarization and the atomic fraction 
inside the target cell were continuously measured [[5], ^TJ. The average values of the target 
polarization for both hydrogen and deuterium data are shown in Tab. EL The luminosity 
was measured by detecting e + e~ pairs originating from Bhabha scattering of the beam 
positrons off electrons in the target atoms, and also 77 pairs from e + e~ annihilations [42]. 



Year 


Target 


Luminosity 


Average Beam 


Average Target 






pb- 1 


Polarization 


Polarization 


1996 


H 


12.6 


0.528 ±0.018 


0.759 ± 0.032 


1997 


H 


37.3 


0.531 ±0.018 


0.851 ± 0.032 


2000 


D 


138.7 


0.533 ±0.010 


0.846 ± 0.030 



Table 1: Integrated luminosities, average beam and target polarizations for the data used in this 
analysis. 



The Hermes spectrometer [43] consisted of two identical halves separated by a hor- 
izontal flux diversion plate, which limited the minimum detected angle. The geometrical 
acceptance was ±170 mrad in the horizontal (bending) plane and between ±(40 — 140) mrad 
in the vertical plane resulting in a range of polar angles between 40 mrad and 220 mrad. 
Each half was instrumented with 3 planes of hodoscopes, 36 planes of drift chambers, 
and 9 planes of proportional chambers. The particle identification system consisted of 
an electromagnetic calorimeter, a pre-shower hodoscope, a transition-radiation detector, 
and a Cerenkov detector. Detailed descriptions of these components can be found in 
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Refs. [43, 44, ^5|, 46, 47, 48]. Positrons within the acceptance could be separated from 
hadrons with an efficiency exceeding 98% and a hadron contamination of less than 1%. 

The main Hermes physics trigger was formed by a coincidence of hits in the hodoscopes 
in the front and back regions of the spectrometer with the requirement of an energy deposit 
above 1.4 GeV in the calorimeter. This trigger was almost 100% efficient for positrons with 
energies above threshold. Events with no positron in the acceptance were recorded using 
a mixture of the main trigger and another trigger formed by a coincidence between the 
hodoscopes and two tracking planes, requiring that there is at least one charged track. The 
influence of trigger efficiencies on the analysis has been studied in [33]. 



3. Experimental Results 

The ratio ^2(x,/i 2 ) of helicity-dependent to helicity-averaged gluon distributions, i.e. the 
gluon polarization, is determined by measuring the double-spin cross section asymmetry of 
one or two high-py inclusive hadrons produced in the scattering of longitudinally polarized 
positrons incident on the longitudinally polarized target. The definitions of the kinematic 
variables in electroproduction used in this paper are shown in Tab. ||. The longitudinal 
double-spin cross section asymmetry is given by the ratio of helicity-dependent to helicity- 

— > — > — y — ► 

averaged cross sections A=Aa/(2a), where a = (a <= + a = *)/2, Aa =a^— a^, and the 
single (double) arrows denote the relative alignment of longitudinal spins of the lepton 
(nucleon) with respect to the lepton beam direction. 

The data for this analysis were collected in 1996, 1997, and 2000 (see Tab. §. The 
analysis presented in this paper includes all (unidentified) charged hadrons. Separate 
asymmetries are given for each charge, target, and event category. 

3.1 Event categories 

Simulations indicate that subprocesses involving hard gluons are relatively enhanced by 
measuring hadrons with high pr with respect to the virtual photon direction {pt{^*))- Cor- 
relations between hadrons in an event may also enhance the signal. Events are categorized 
by the number of hadrons observed in an event and whether kinematic information on the 
scattered positron is available or not. Each possible combination of two hadrons is counted 
as a separate event in the pairs category. The categories are defined in detail as follows: 

• 'anti-tagged' single inclusive hadrons: Events with leptons in the acceptance 
were not included in this category. The hadron transverse momentum PT(beam) was 
measured with respect to the beam direction as the direction of the virtual photon is 
unknown. In most cases, the undetected positron had a small scattering angle (and 
hence Q 2 is small) and stayed inside the beam pipe. The difference between PTtbeam) 
and Pt(-j*) is then very small. However, the positron could also escape the detector 
acceptance because of a large scattering angle, in which case Q 2 was large. The large 
angle of the virtual photon with respect to the beam axis results in a significantly 
larger PT(beam) than Pt(-j*) of the hadron. Although these events with large Q 2 are 
rare, they can account for a significant fraction of the hadrons at high PT(beam)- F° r 
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k = (E, k), k! = (£", k') 4-momenta of the initial and final state leptons 

9, (ft Polar and azimuthal angle of the scattered positron 

M Mass of the initial target nucleon 

q = (E — E' ,k — k') 4-momentum of the virtual photon 

Q 2 = —q 2 AEE' sin 2 | Negative squared 4-momentum transfer 

v = ^jj- l = E — E' Energy of the virtual photon 

x = 2T^g ^ 1 HIu Bjorken scaling variable 

y = ^| l = -| Fractional energy of the virtual photon 

W 2 = {P + qf 

l = M 2 + 2Mu — Q 2 Squared invariant mass of the virtual-photon nucleon system 

p = (Ehjp) 4-momentum of a hadron in the final state 

Pt Transverse momentum of a hadron 

Pt(j*) Pt with respect to the virtual photon 

PT(beam) Pt with respect to the incoming positron 

Transverse hadron momentum from fragmentation 

£ P 2 T(beam) For two hadrons: (j)f {beam) ) 2 + {pf [beam) ? 

z = ^p^j- l = ^ Fractional energy of the final state hadron 

x Parton momentum fraction 

s = (p a + pi,) 2 Mandelstam variable for partonic process ab — > cd 

t = {p a — p c ) 2 Mandelstam variable for partonic process ab — > cd 

u = (pb — Pc) 2 Mandelstam variable for partonic process ab — s- cd 

pQCD scale 

Pt (= \l ^ for m = 0) Transverse momentum of final state partons 

in the CM-system of the hard subprocess 

kx Intrinsic transverse momentum of partons 

in the nucleon and photon 

Table 2: Definition of kinematic variables. 



Pt > 1-0 GeV the deuteron (proton) data sample in this category contains 1272k 
(419k) hadrons. 

• 'tagged' single inclusive hadrons: The scattered positron has been detected with 
Q 2 > 0.1 GeV 2 , W 2 > 4 GeV 2 , and y < 0.95. The hadron transverse momentum 
Pt(j*) is measured with respect to the virtual photon direction. For pt > 1 GeV this 
deuteron (proton) data sample contains 53k (19k) hadrons. 

• inclusive pairs of hadrons: The hadron pair sample consists of all pairs of charged 
hadrons with PT(beam) > 0-5 GeV. The transverse momentum PT(beam) is measured 
with respect to the beam direction, because only in 10% of the events the positron 
was detected. With the additional requirement X]PT(beam) > 2.0 GeV 2 the deuteron 
(proton) data sample contains 60k (20k) hadron pairs. With this requirement applied, 
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6% of the anti-tagged inclusive hadrons with PT(beam) > 1-0 GeV are contained within 
the pairs sample. 

For all three categories, Hermes data are available for various combinations of target 
and/or hadron charge detected. As the samples differ in the hard subprocess and final 
state kinematics and fractions of contributing subprocesses, the corresponding results for 
the gluon polarization -^-(x, fJ- 2 ) provide a measure of the consistency of the extraction. The 
final result for ^-(x, fi 2 ) is obtained from the anti-tagged inclusive hadrons originating from 
a deuterium target. The other data samples have too small a statistical power to justify 
carrying out the extensive analysis needed to obtain the systematic uncertainties. 

3.2 Asymmetry results 

The double-spin asymmetry measured is given by 

A-meas = All — — — — — . (3-1) 

— > — ► 

Here (N^) is the number of hadrons or hadron pairs for target spin orientation parallel 

— > — > 

(anti-parallel) to the beam spin orientation, Lr^ \L ) is the corresponding integrated 

— > — ► 

luminosity, and Up (Lp) is the integrated luminosity weighted by the live-time fraction 
and the absolute values of beam and target polarizations. There is a small background 
(< 0.1%) from positrons misidentified as hadrons (and vice versa). In the tagged category 
a correction was applied for an approximately 5% contribution of positrons originating 
from charge-symmetric processes. 

The asymmetries for the anti-tagged and tagged categories are shown as a function of 
transverse momentum in Figs. || and |3|, respectively, and listed in tables ^ - |l0[ The asym- 
metry of the pairs is presented as a function of the minimum requirement, (^ PT(beam)^ min > 
in Fig. |] and in table [ll]. The considerably different values of the asymmetries in the 
different categories, charges and targets are due to the different underlying mixtures of 



subprocesses and of quark content, as discussed in Sect, p.5 . 

The curves in Figs. ||, ||[ and |I| show the asymmetries calculated by the procedure 
discussed in Sect. ||] using the values ^(x,// 2 ) = —1, 0, +1 (from top to bottom). They 
illustrate the sensitivity of the Hermes data to ^f-( x , f^ 2 )- The data are close to the central 
curve indicating small average values of 

4. Physics Model 
4.1 Subprocesses 

Both the helicity-averaged and helicity-dependent cross sections include contributions from 
hard subprocesses that can be calculated using pQCD and from soft subprocesses such as 
those described by the Vector-Meson Dominance (VMD) model (see Fig. [l]). A smooth 
transition from soft subprocesses to hard subprocesses is regulated by a set of cutoff pa- 
rameters (for details [33, |4^, ^]). The measured asymmetry is the weighted sum of the 
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Figure 2: Measured asymmetry for the anti-tagged category of events for positive (left) and 
negative (right) inclusive hadrons from hydrogen (top) and deuterium (bottom) targets as a function 
of PT(beam)- The uncertainties are statistical only. There is an overall normalization uncertainty of 
5.2% (3.9%) for hydrogen (deuterium). The curves show the Monte Carlo asymmetries for three 
different fixed values assumed for the gluon polarization. 



asymmetries of all subprocesses. When it is impossible to reliably separate the subprocesses 
experimentally, as in fixed-target experiments, the fractions of events originating from the 
different subprocesses must be modeled. In the analysis described in this paper, this is 
done using the spin independent Monte Carlo (MC) program Pythia 6.2 [49, 5C, 51]. 

The various subprocesses are classified in terms of the model used in Pythia. In 
this model, the wave function of the incoming photon has three components, a "VMD", 
a "direct" and an "anomalous" one. The generic photon processes following from this 
decomposition are depicted in Fig. [|. The direct photon interacts as a point-like particle 
with the partons of the nucleon, while the VMD and anomalous components interact 
through their hadronic structure. 

Figure |5]b shows an example of a direct process. The direct pQCD subprocesses studied 
in this analysis are the O(o®) DIS process (Fig. [j]a), the 0(a).) processes PGF (Fig. |l]b), 
and QCDC (Fig. [TJc) . 

The VMD component is characterized by small-scale, non-perturbative fluctuations 
of the photon into a qq pair existing long enough to evolve into a hadronic state before 
the interaction with the nucleon. This process can be described in the framework of the 
VMD model, where the hadronic state is treated as a vector meson (e.g., p°, u, <p) with the 
same quantum numbers as the photon. Higher-mass and non-resonant states are added in 
the Generalized VMD (GVMD) model. The (G)VMD hadronic states can undergo all the 



-9- 



0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 



L ep->eh + X . 
» »■«-.-,-.-» — I — ; 




L ep->eh~X-' 


»»* 




> 


i m, A i 




r MC, Ag/g= -1 

r MC, Ag/g= 

L MC, Ag/g= +1 


, 


• "■•v;-l ; 

-_ \ 

L • Hermes 

: 1 i , , , 1 , , , , 1 , , , , 1 , 


1 

. " .- -J . . 


L ed^eh + X 

— *> "* i 

: !».r».f ftoaC * * * . . . .73 


— — 

(I 

I -. 


L ed->eh"X , 

-_ / 
— * 


1 > 


~. ... \ .... \ \ , . 




~\ .... \ .... \ .... \ . 


1 

, , 1 , 



0.5 1 1.5 2 0.5 1 1.5 2 

P T(f) [GeV] 

Figure 3: Measured asymmetry for the tagged category of events for positive (left) and negative 
(right) inclusive hadrons from hydrogen (top) and deuterium (bottom) targets as a function of 
Pt(7*)- The uncertainties are statistical only. There is an overall experimental normalization un- 
certainty of 5.2% (3.9%) for hydrogen (deuterium). The curves show the Monte Carlo asymmetries 
for three different fixed values assumed for the gluon polarization. 
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Figure 4: Measured asymmetry for hadron pairs produced from hydrogen (left) and deuterium 
(right) targets as a function of the minimum value of PT(beam) • The uncertainties are statistical 
only. There is an overall experimental normalization uncertainty of 5.2% (3.9%) for hydrogen 
(deuterium) . The curves show the Monte Carlo asymmetries for three different values assumed for 
the gluon polarization. 



interactions with the nucleon allowed in hadronic physics, i.e., elastic and diffractive as well 
as inelastic non-diffractive reactions. The latter can be either soft ("low-pr") processes or 
hard QCD 2 — > 2 processes. A generic example of a VMD process is shown in Fig. [5]a. 
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a) b) c) 

Figure 5: Contributions to hard p interactions: (a) VMD, (b) direct, and (c) anomalous. Only 
the basic graphs are illustrated; additional partonic activity is allowed in all three processes. The 
presence of spectator jets has been indicated by dashed lines, while full lines show partons that 
(may) give rise to high-pt jets. The gray ovals represent multiparton wave functions. Anomalous 
states are built up from a perturbatively given qq fluctuation, while VMD fluctuations allow no 
simple perturbative representation; hence the difference in the placing of the ovals [BU . 



The anomalous photon is characterized by sufficiently large-scale, perturbative fluc- 
tuations of the photon into a qq pair. The allowed processes are the same pQCD 2 — >■ 2 
processes as in the hard VMD case, with the difference that for the anomalous component 
the parton distributions of the photon are relevant, whereas for the description of the hard 
VMD component those of the vector meson are used. Both hard VMD and anomalous 
components are usually referred to as "resolved" photons. Depending on whether a quark 
or a gluon in the nucleon is struck by a resolved photon the corresponding hard subpro- 
cesses are labeled by a 'q' or 'g' in this paper. A generic example of an anomalous process is 
shown in Fig. ||c. The resolved-photon processes are of 0(a).) with a hidden l/a s (fi 2 ) term 
in the evolution of the photon's parton distributions canceling the additional vertex |52]. 

For hard subprocesses the nucleon is described by helicity-averaged (helicity-dependent) 
PDFs, which are the average (difference) of the number densities of partons of type f whose 
spins are aligned, / + , those whose spins are anti-aligned, /~, with respect to the nucleon 
spin: f(x,n 2 )=f+(x,fi 2 ) + f'(x,^i 2 ) (Af(x,^)=f + (x,fi 2 ) - f-(x,fi 2 )), where /= u, d, 
s, or g. The integral over x, A/(/i 2 ) = tidxAf(x,[i 2 ), gives the total spin contribution 
of the respective partons to the nucleon spin, as used in Eq. |Ll]. The hard part of the 
single-inclusive differential helicity-dependent cross section for the process ^y*p — > hX can 
be expressed as an integral over the parton distribution functions, the hard partonic cross 
sections for the subprocesses ab — > cX, and the fragmentation functions. It can be written 
schematically as 

dAa r P ^hX = J2 J dx a dx b dz c Arf(x a ,^ 2 )Af b N (x b ,fi 2 ) 

a,b,c=q,q,g 

KdAa ab ^ x {sXv\Q 2 )D h c {z c ,ii 2 ) , (4.1) 

and correspondingly for the helicity-averaged cross section and distributions. Here x b is the 
fraction of the nucleon momentum carried by parton b and f b (x b , fi 2 ) (A/^ v (xf ) , ^ 2 )) is the 
corresponding nucleon PDF. Similarly x a is the fraction of the photon momentum carried by 
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parton a, and f2 (x a ,fi 2 ) (A/a (x a ,/x 2 )) is the corresponding photon PDF. For the direct- 
photon processes a equals 7* and f2 (x a ,n 2 ) (A/J (x a ,fi 2 )) reduces to 5(1 — x a ). The 
fragmentation function D^(z c ,fi 2 ) describes the hadronization of a parton c into a hadron 
h with a momentum p^ = z c p c . The hard partonic cross section da ab ^ cX (s,i, fJ?,Q 2 ) 
(dAa a ~* cX (s, i, /i 2 , Q 2 )) depends on the subprocess kinematics, the renormalization and 
factorization scales, and on Q 2 in case of the direct-photon processes. Here, s and t are the 
Mandelstam variables for the partonic interaction, which are related to x a and Xb- More 
information on the kinematic variables is given in table 0. The cross section for hadron 



pairs da~ f * p -> hlh2X (dAa~<* p -> hlh2X ) can be obtained analogously to Eq. [O. 

The cross sections and asymmetries of the soft VMD interactions can only be mod- 
eled phenomenologically. The Pythia model incorporates the total 7^ and hadronic cross 



section parameterizations of Donnachie and Landshoff |53j together with quark counting 



rules [54, 55]. This model successfully describes the measured total, elastic, and diffractive 
cross sections over a wide energy range. The non-diffractive cross section is modeled in 
Pythia as the difference of the total cross section and the summed elastic and diffractive 
cross sections; the corresponding subprocess is called "low-p^". The Pythia model pro- 
vides a smooth transition from real to virtual photons and is applicable from very small to 
large values of Q 2 . It uses a number of cutoff, scale, and suppression parameters together 
with several possible prescriptions on how to use them to select the underlying subprocess 
of an event. The default prescriptions and the cutoff and scale parameters were developed 
and tuned to match high energy data. In this application to the lower energy of Hermes 
the influence of various prescriptions and parameter values has been carefully studied (see 
Sects. |and 0). 

Table |3| shows a compilation of the modeled reactions, the corresponding Pythia 
subprocess numbers, their classification, description, and name used in this paper. 

4.2 Signal and Background Asymmetries 

In the simulation, the cross section is considered to arise from an incoherent superposition 
of all contributing subprocess amplitudes. The kinematic selection criteria (e.g., event cate- 
gory and hadron pr) for the Monte Carlo are the same as those for the data. Pythia events 
are generated independent of helicity, therefore the MC asymmetry Amc is calculated by 
weighting each selected MC generated hadron with the calculated event asymmetry Wk- 
The average of these weights is Ai, the asymmetry for subprocess i 

1 Ni 

^ = n:J2 w ^ ( 4 - 2 ) 

* fe=i 

where Ni is the number of entries. The event-by-event weighting method guarantees the 
correct integration over the subprocess kinematics, and all partons in the nucleon and 
the photon (where applicable). The Monte Carlo asymmetry Amc is the sum of the 
asymmetries from signal (Afj G ) and background (Af G c ) subprocesses weighted by their 
fraction of entries R SIG and R BG . It is given by 

A MC (pt) = R BG Af G ( PT ) + R sig AHS(pt) = £ RiA + Ri^i, (4.3) 

i£BG iGsig 
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Subprocess 


# 


Class 


Description 


Name 


soft VMD 


VN 


-»■ FiV 


91 


background 


elastic VMD 


exclusive VMD 


VN 


-»■ FX 


92 


background 


single-diffractive VMD 




VN 


-»■ XX 


93 


background 


single-diffractive VMD 




VN 


-> XX 


94 


background 


double-diffractive VMD 




VN 


-> X 


95 


background 


soft non-diffractive VMD 


\ow-px 


RESOLVED (hard VMD and anomalous) 


qq -> 


99 


11 


background 


QCD 2^2 


QCD 2 -> 2(g) 


99 -> 


99 


12 


background 






99 -> 


99 


13 


background 






59 -> 


59 


28 


background 






95 -> 


95 


28 


signal 




QCD 2 -> 2(g) 


55 


99 


53 


signal 






55 -> 


55 


68 


signal 






DIRECT 


* 

7 9 - 


■+ 9 


99 


background 


LO DIS 


DIS 


7t9" 


-> 95 


131 


background 


(transverse) QCDC 


QCDC 


* 

7l9" 


->■ 95 


132 


background 


(longitudinal) QCDC 




* 

7t5 - 


-> 99 


135 


signal 


(transverse) PGF 


PGF 


* 

7l5 - 


-> 99 


136 


signal 


(longitudinal) PGF 





Table 3: Description of the subprocesses used in this paper. Columns from left to right: subprocess, 
Pythia subprocess number, classification as signal or background, description, and name used in 
this paper. A vector meson is denoted by V. 



where Ri is the fraction of entries from the subprocess i calculated in the PYTHIA simu- 
lation. Background processes are all subprocesses that do not involve a hard gluon from 
the initial nucleon. These include all soft processes, the direct processes DIS and QCDC, 
and all resolved pQCD processes, which involve a quark or antiquark in the nucleon, i.e., 
QCD 2 — > 2(g). They are listed in Tab. ||. All subprocesses involving a hard gluon of the 
nucleon in the initial state are considered to be signal processes, i.e., PGF and the hard 
2 — > 2(g) processes. 

The event-by-event weight uu for hard subprocesses is given by 



A/f ( w 2 ) Af b N (x b ,» 2 ) 
fT(x a ,^) ' f b N (x b ,^ 



w = a{s,t,» ,Q )• . — • N( 2 . , (4.4) 



where A /J /fa = 1 for x a = 1, i.e., direct photon processes. The hard subprocess 
asymmetry is a(s,t, /i 2 ,Q 2 ) = Aa/(2a). The lowest order equations for important hard 
subprocess asymmetries are compiled in appendix [A]. The VMD and GVMD diffractive 



subprocesses may have small asymmetries at Hermes energies [gg, 57, 58, |5S|]. The asym- 
metry of the 1ow-pt process was estimated from the measured asymmetries and found to 
be non-zero (see Sect. 5.3). In both cases, the virtual-photon depolarization factorD(y, Q 2 ) 
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(see Eq. A.15| ) has to be applied to the weight in order to account for the transformation of 



the virtual-photon nucleon asymmetry into a lepton-nucleon asymmetry. The asymmetry 
from signal subprocesses depends on the unknown ^-(x,p 2 ) averaged over the subprocess 
kinematics in the specified px range. It can be written as 

N S1G 



N SIG 

k=l 



Af2\x a ,fi 2 ) Ag 2 \ SI " 
— — — • —(x b ,n )) 

fa [X a ,p 2 ) 9 / 



a(s,t,n ,Q ) • — 3T- — (x b ,fj, ) ) (pr), (4.5) 



where N SIG is the number of entries from all signal processes. The extraction of the 
quantity of h 
A mc(pt) by 



quantity of interest, ^2(x,/U 2 ), is based on Eq. L5 replacing the unknown asymmetry 



\SIG ( ^ _ A meas ( PT ) ~ R BG Aff c (p T ) 



A" u (pr) = m ^ ± ' RSIG . (4.6) 



In Sect. || methods will be described to extract ^ from the right hand side of Eq. 4.5 



5. Monte Carlo simulation 

The relevant subprocess cross sections have been modeled by the Pythia Monte Carlo 



program, which uses Jetset [6C] for describing the fragmentation process. The standard 



helicity- averaged input PDFs used are CTEQ5L [61] for the nucleon and Schuler and 
Sjostrand |62| for the photon. The scale p? of the 2 — >• 2 subprocesses is defined to be 



fi 2 = p\ + \Q 2 (also commonly referred to as Q 2 ). Electromagnetic radiative effects [63, 64] 
have been added to Pythia and they constitute a relatively small correction for hadron 



production at Hermes kinematics [33]. Events generated by Pythia are passed through 



a complete Geant 3 ]63] simulation of the Hermes spectrometer. 
5.1 Monte Carlo tuning 

In order to account for the relatively low center-of-mass energy of the Hermes experiment 
several parameters in the event generation were adjusted and the model describing exclusive 
vector meson production was improved 33]. This was done in the kinematic region of the 



tagged events because more kinematic variables are measured for this category than for 
the anti-tagged category. The tuning of the fragmentation parameters [36] was performed 



using a subsample with PT(-y*) < 0.8 GeV and Q 2 > 1 GeV 2 where the DIS process (Fig. |l|a) 
is dominant and NLO corrections are small. The values of the adjusted parameters, shown 
in Tab. |7| in appendix [l|, are used for all event categories. 

Figure ^ shows the measured and the simulated cross sections as a function of xb, Q 2 , 
and z for the tagged category of events using a deuterium target. Both the simulated and 
measured cross sections are not corrected for acceptance effects. These cross sections vary 
over more than three orders of magnitude. The data and MC simulation agree to within 
15% for xb < 0.2, where most of the data reside for the tagged event category. Thus in 
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Figure 6: Top panels: Measured cross section in the Hermes acceptance for tagged hadrons as 
a function of xb (left), Q 2 (middle), and z (right) for positive (full points) and negative hadrons 
(open points) using a deuterium target. The lines show the tuned Pythia 6.2 calculation. Bottom 
panels: The corresponding ratios of the Pythia calculation to the measured cross section. 



this region the modified Pythia 6.2 program with the adjusted parameters gives a good 
representation of the cross section at Hermes energies. 

The description of the kinematic dependences of the tuned Monte Carlo code for the 
individual subprocesses must be consistent with independent LO pQCD calculations |67| . 
Such calculations presently exist only for inclusive 7T° production and only in the collinear 
approach, where the intrinsic transverse momentum /ct of the partons in the nucleon and in 
the virtual photon, and also the transverse momentum p^ ag arising from the fragmentation 
process are set to zero. 

For a comparison of Pythia with these LO pQCD calculations a special simulation 
with kx = and p^ 0-9 = was performed, by replacing the string fragmentation performed 
by JETSET with weights obtained from the fragmentation functions of Ref. j68[| . The 
resulting transverse momentum PT(beam) of the 7T° is calculated according to PT(beam) = 
z ■ pt- Both this simulation and the pQCD calculation are performed in the Hermes 
kinematics for inclusive 7r° production at Q 2 < 0.01 GeV 2 , 0.2 < y < 0.9, disregarding the 
detector acceptance. Figure |7] compares the resulting cross sections for resolved photon, 
QCDC, and PGF processes from the simulation and the pQCD calculation. In the collinear 
approach the DIS subprocess is not included, because the Pt(-j*) oi the final state hadron 
is zero, and also for low Q 2 (Q 2 < 0.01 GeV 2 ) it does not result in a sizable PTibeam) ■ 

The agreement between the simulated cross sections for the individual subprocesses 
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Figure 7: Cross sections for inclusive tt° production from resolved photon, QCDC, and PGF 
processes simulated using Pythia (solid lines) compared to the LO pQCD calculations from Ref. |6£j 
(full points). Simulation and calculation are done in the collinear approach at Q 2 < 0.01 GeV 2 , 
0.2 < y < 0.9. Green/Grey lines: subprocess cross sections after varying the renormalization and 
factorization scales by factors of k and 2 in the simulation. 



and the calculations is well within the scale uncertainty (i/x , 2/i ) of the simulation (the 
dashed lines in Fig. |7|). The LO pQCD calculations show a similar dependence on the 
variation of the renormalization and factorization scales {\pt (beam) i^Pt (beam)) i see Fig- H 



in Ref. §7]. 



5.2 Effects of intrinsic and fragmentation transverse momenta 

While the effect of intrinsic and fragmentation transverse momenta cannot yet be studied 
in LO pQCD calculations, a Pythia simulation can be used. For the standard simulations 
presented in this analysis a Gaussian distribution with a 0.4 GeV width is used for both 
hr and p^ a9 |^6|. These values are consistent with those obtained in Ref. [fT0[| . Both 
intrinsic and fragmentation transverse momenta alter the relationships of p^ to PT(beam) > 
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from PT(beam) = z ' Pt to PT(beam) = z (^t + Pt) + 9 > and hence the distribution of 
and x. This in turn influences the dependence of the cross section on PT(beam) ■ The effects 
on the cross section for inclusive 7r° production from the PGF subprocess, of first adding 
nonzero kx = 0.4 GeV and secondly using Jetset with p T rag = 0.4 GeV are shown in 
Fig. |8|. Including only kx in the simulation decreases (p T ) from 1.9 GeV 2 to 1.6 GeV 2 and 
(x) from 0.32 to 0.28, and increases the cross section by a factor of two. Including both 
kx and p T rag further decreases (p T ) to 1.1 GeV 2 and (x) to 0.22, and increases the cross 
section by another factor of 10. These studies show that at fixed-target kinematics, like at 
Hermes, intrinsic and fragmentation transverse momenta cannot be neglected in pQCDC 
calculations. Similar conclusions were drawn in Ref. [fn[]. Perhaps resummation techniques 
1 72], which account for initial and final state radiation effects, can help to achieve more 
realistic calculations. 




Figure 8: The simulated cross section for inclusive ir° production from the PGF subprocess vs. 
PT(beam)i Pj" an d x f° r Q 2 < 0-01 GeV 2 , 0.2 < y < 0.9. The simulations are done using the 
collincar approach ((green/gray) solid line), the collinear approach together with intrinsic kx for 
the partons in nucleon and photon (dashed line) , and with intrinsic kx together with fragmentation 
transverse momenta (solid line). For the first two simulations fragmentation is modeled using the 
KKP-fragmentation functions f6q| , for the third one Jetset with the standard settings listed in 
Table is used. 



5.3 Analysis of Monte Carlo events 

Pythia events are used to calculate cross sections, individual subprocess fractions Ri 
and event weights (w)i within the Hermes acceptance. The event weights for the pQCD 
processes (Eq. (P[ ) are obtained using the hard subprocess asymmetries (see appendix 
|A|) and GRSV (standard scenario) ]73| helicity-dependent PDFs in conjunction with the 
GRV98 Jrlj helicity-averaged PDFs to calculate Af N /f N for the nucleon. In order to 
calculate A/ 7 // 7 for the photon the averages of the maximal and minimal scenarios of the 
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GRS [f7f| |76|] helicity-dependent PDFs are used in conjunction with the GRS [77] helicity- 
averaged PDFs. 



For elastic and diffractive VMD processes the asymmetry is set to zero [59|. For the 
low-px process two alternative assumptions for the asymmetry have been investigated: 
Ai ow - PT = D(y, Q 2 ) ■ A\ and A 

low-pr — where A± is a parameterization of the photon- 
nucleon asymmetry in inclusive DIS. The resulting MC asymmetries and the measured 
asymmetry are shown in Fig. ^ (left) vs. xb for the tagged category and a hydrogen 
target. The corresponding deuterium data are not shown because for this target both 
assumptions are indistinguishable and match the data. For the anti-tagged category (Fig. || 
(right)) the PT(beam) dependence of the double-spin inclusive asymmetry Aii is shown for 
both targets. The model A\ ow - PT = D ■ A\ matches the data better than Ai ow - Pt = in the 
kinematic domains where R 1ow 'Pt \ s i ar g e (\ow-xb for tagged and the lowest px for anti- 
tagged categories, respectively) and contributions of hard QCD processes are negligible. 



The standard asymmetry for the low-pr process was chosen to be Ai ow . 



PT 



D-At, because 



of this agreement and because the semi-inclusive asymmetry for all charged hadrons is 
approximately equal to the measured inclusive asymmetry. The world data on A\ have 
been parameterized by a + x b B ■ (1 — e CXB ) for xb > 10~ 3 , and extrapolated to the smaller 
xs-values {{xb) ~ 0.0001) typical for the anti-tagged sample. 
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Figure 9: Left panels: The measured double-spin inclusive asymmetry A\\ (full points) and two 
MC asymmetries (solid and dashed lines) based on different assumptions for the low-pr subprocess 
asymmetry (top) and the process fractions vs. xb (bottom) for tagged events on a hydrogen target. 
Right panels: Double-spin asymmetry vs. PT(beam) for anti-tagged hadrons on a hydrogen (top) 
and deuterium (bottom) target. 



To avoid any bias from the experimental trigger to the results presented, the MC 
events received an additional weight to account for trigger inefficiencies, if measured and 
simulated cross sections are compared. The pt dependences of the cross sections, individual 
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subprocess fractions Ri, average event weights A{, and weighted asymmetries R%Ai for the 
three event categories are shown in Figs. [l0| (anti-tagged), |ll| (tagged), and |l2] (pairs). 
All three categories have in common that: 

• The cross sections span four orders of magnitude, decreasing rapidly with pt\ 

• Reasonable agreement between data and Monte Carlo is observed for low transverse 
momenta. With increasing pt the Monte Carlo description becomes worse, underes- 
timating the data by up to a factor of four at the largest pt; 

• The fractions R 1ow ~pt anc j pexcl.VMD decrease with increasing px and the corre- 
sponding asymmetries are very small or zero, respectively; 

• In general the contributions from hard QCD subprocesses increase with increasing 
PT- At high pt subprocesses involving quarks in the nucleon contribute less than the 
signal processes; 

• The asymmetries for the two signal subprocesses, QCD 2 — > 2(g) and PGF have 
opposite sign. For a positive gluon polarization like that of GRSV, this results in 
a sizable negative asymmetry for PGF, and positive asymmetries for the 2 — > 2(g) 
processes; 

• Some asymmetries and fractions depend on the charge of the hadron. 

Even though soft effects from initial and final state radiation and additional nonper- 
turbative processes are taken into account in the Pythia simulation, the Monte Carlo 
simulation still fails to describe the cross sections at pt > 1 GeV. This shortcoming may 
be explained by missing large higher order corrections to the hard processes. These correc- 
tions have been evaluated for the next to leading order (NLO) cross section in Ref. ||S7f| , in 
the collinear approach for Q 2 < 0.01 GeV 2 and pj? 1 = zpT > 1 GeV, for all hard processes 
(QCD 2 — > 2, PGF, QCDC) contributing in this region. The similar kinematics of hard 
processes in the pQCD-calculation and the PYTHIA simulation allows one to approximate 
the effect of NLO corrections to the Monte Carlo cross section. A /c-factor, i.e., the ratio 
of LO to NLO cross sections is applied as a weight to each hadron originating from a hard 
process. The /c-factors from Ref. |37| are very large (almost 5) at pt ~ 1 GeV and decrease 
with pt to about 2.5 at pt = 2.4 GeV. For the reweighting of the Monte Carlo events they 
have been extrapolated down to PT(beam) = 0.8 GeV, and it was assumed that p^ 1 can be 
approximated by PT(beam) ( see the discussion in Sect. [O] about the collinear approxima- 
tion). The results shown in the cross section ratio of Fig. |l(] indicate that the inclusion 
of NLO effects to the Monte Carlo could significantly improve the description of the cross 
section. Effects of similar size may exist for the other categories, but NLO calculations for 
those are not yet available. The /c-factors for the asymmetry have also been calculated in 



1 67] and are approximately 2 in the experimental range. Unfortunately it is not possible to 



consistently take into account /c-factors in the extraction of ^-(x,fi 2 ), therefore the result 
will essentially be a LO result subject to potentially large NLO corrections. 
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Figure 10: The cross sections and subprocess contributions, in the Hermes acceptance as a 
function of PTfbeam) f° r the anti-tagged category of events and a deuterium target (left: positively 
charged hadrons, right: negatively charged hadrons). Top: The measured cross section and that 
generated by Pythia. Second row: The ratio of these two cross sections. Also shown is the effect 
of the fc-factor based on Ref. |37j (see text). Third row: The subprocess fractions from Pythia. 
Bottom two rows: The asymmetries and the asymmetries weighted with the subprocess fractions 
for each subprocess using Refs. |73] and []74) for the gluon PDFs. 
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Figure 11: The cross sections and subprocess contributions, in the Hermes acceptance, as a 
function of Pt( 7 *) for the tagged category of events and a deuterium target (left: positively charged 
hadrons, right: negatively charged hadrons). Top: The measured cross section and that generated 
by Pythia. Second row: The ratio of these two cross sections. Third row: The subprocess 
fractions from Pythia. Bottom two rows: The asymmetries and the asymmetries weighted with 
the subprocess fractions for each subprocess using Refs. ]73| and j74j for the gluon PDFs. 
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Figure 12: The cross sections and subprocess contributions, in the Hermes acceptance, as a 
function of the minimum value of 2^,PT(beam) f° r ^ ne production of inclusive hadron pairs on a 
deuterium target. Top: The measured cross section and that generated by Pythia. Second row: 
The ratio of these two cross sections. Third row: The subprocess fractions from Pythia. Bottom 
two rows: The asymmetries and the asymmetries weighted with the subprocess fractions for each 
subprocess using Refs. [JT3] and j74j for the gluon PDFs. 
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For the anti-tagged category the LO DIS fraction dominates the yield of positive 
hadrons at high PT(beam) ■ This is due to the subsample of events with the positron having 
a large scattering angle and missing the Hermes acceptance. The subprocess fractions for 
LO DIS and QCDC are larger for positive hadrons because of M-quark dominance. Both 
signal subprocesses contribute approximately 20% to the cross section at high PT(beam) ■ The 
pairs category has a larger signal fraction than the other categories, but a much smaller 
number of events. The mixture of the background processes and their contribution to the 
background asymmetry is different for each category. 



6. Determination of the gluon polarization 



6.1 Kinematic considerations and requirements 

The average value of in a pt range is determined directly from Eq. |6.l| (see Sect. 



However, as shown in Fig. 13J, there is a large range of x spanned by the data for each pt 
range. In order to circumvent this difficulty, the value of ^ and the appropriate value of 
x is determined through a minimization procedure using a functional form for ^(x) (see 
Sect. |6.3| ). The scale dependence of ^ is neglected because almost all pQCD models are 
monotonic and vary slowly as a function of fj? over the relatively small relevant range. In 
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Figure 13: The range of generated x for different values of PT(beam) calculated by Pythia for all 
signal processes, for the anti-tagged category of events and a deuterium target. 



order to optimize the accuracy of — the following criteria that maximize the sensitivity 
of the MC asymmetry to are applied to the individual data samples: 
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These requirements balance the statistical accuracy of the measured asymmetries (decreas- 
ing with Pt-, as shown in Figs. §-||) against the signal process fractions (increasing with p?, 



as shown in Figs. 10-12). For the events within these limits it is observed that: 
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• The Pythia simulations displayed in Fig. |TJ show a strong correlation between the 
hard scattering transverse momentum (p T ) of the signal subprocesses and the mea- 
sured hadronic px (%2Pt)\ 

• For larger values of pr, there is greater sensitivity to the hard processes involving 
the gluon (see Figs. |2||4|), which leads to reduced systematic uncertainties due to 
corrections for background asymmetries. 
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Figure 14: The correlation of the average hard scattering (p\) of all signal subprocesses as 
calculated by Pythia with the hadron pr for inclusive hadrons as calculated for the experimental 
data for the deuterium target. Left: tagged category; Center: anti-tagged category; Right: hadron 
pairs category. The dotted line goes along (p^) = pj*(h) ((Pt) — J2Pt/^) an d the vertical dashed 
line shows the minimum px (%2Pt) use d for the analysis. 



The gluon polarization — is determined using Eqs. 4.3, 4.5, and 4.6. The anti-tagged 



category has sufficient statistics to allow extraction of ^ in four PT(beam) bins (1.0 - 1.2 - 
1.5 - 1.8 - 2.5 GeV), which are obtained by combining the bins shown in Fig. |2|and table |9|. 
The other categories are represented by a single range in px- 

6.2 pt dependence of ^ 

If the dependence of — (x,/u 2 ) on x and /x 2 is weak in the limited kinematic range of the 



experiment, can be factored from the r.h.s of Eqs. 4.5, so that together with Eqs. 

4.6 we obtain for the gluon polarization averaged over the covered x and /j? ranges 
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(6.1) 



where the subprocess fractions and kinematics are determined using Pythia. As is shown 
in Fig. [l^, different ranges in pjn correspond to different ranges and distributions in x. It is 
intrinsic to this method that there is no knowledge on the dependence of on x, therefore 
no meaningful value of the average x can be determined by this method, which nevertheless 
can be used as a consistency check between the different independent data sets. 
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The results for different event categories, targets and hadron charges are listed in 
table |] and shown in Fig. |l5] as a function of pr- The results for the pairs category 
are displayed at the average y ^p|,/2, and those for the tagged category at the average 
Pt(~{*)- Each of these data sets has a somewhat different mixture of background and signal 
processes as a function of pr, as seen in Figs. The measured values of (^)(pr) 

should be equal for both targets and both hadron charges because of the same range in 
x and /i 2 . The values shown in Fig. 15 within each category and for each PT(beam) bin 
indeed agree in general within the statistical uncertainties. This is a strong indication that 
Pythia provides a consistent description of the underlying physics. The systematic charge 
dependence is accounted for by assigning a systematic uncertainty to the value of the p T mfil 
(Pythia parameter PARJ(21)). 

In |2^] the kinematic selections for the hadron pairs used to calculate the asymme- 
try to extract ^0*0 was > 1.5 GeV and p T 2 > 0.8 GeV. These events are mostly 
contained in the event sample used to calculate the asymmetry in the left panel of fig. [| 
if £ (P 2 T(beam)) min = P^eam)? + ^T%e am )) 2 > 3 -° GeV2 is required. The asymmetry 
for hadron pairs with V (pL p „„>) > 3.0 GeV 2 presented here is statistically consistent 

with the average asymmetry for p^ 1 > 1.5 GeV and p T 2 > 0.8 GeV from p9| . The difference 
for (^■}(pt) obtained for the inclusive pairs of hadrons in this paper compared to the result 



presented in [£9| can be explained by the different treatment of the underlying signal and 
background subprocesses contributing to the asymmetry and the difference in kinematic 
selections of the hadron pairs used calculating the asymmetry. The model presented in 
1 29] used only 2 subprocesses (PGF and QCDC) to describe the measured negative asym- 
metry for the proton target. For the determination of the subprocess fractions also the 
VMD process was considered, which was treated to have no subprocess asymmetry, which 
is consistent with the model used in this paper. The resulting subprocess fraction for PGF 
in [^] is bigger than from the model presented in this paper. This combined with the pos- 
itive asymmetry for the QCDC subprocess leads to the sizable positive gluon polarization 
reported in p9| (Note: a(s,t, fi 2 ,Q 2 ) is negative for PGF in the probed kinematics). 





Proton 


Deuteron 


(pr) (GeV) 


h+ 


h~ 


h+ 


h" 


anti-tagged 


1.11 


-0.076 ±0.150 


0.201 ±0.162 


-0.063 ± 0.096 


0.125 ±0.096 


1.30 


0.011 ±0.120 


0.125 ±0.103 


-0.005 ± 0.073 


0.080 ± 0.059 


1.60 


0.116 ±0.195 


0.619 ±0.174 


-0.087 ±0.119 


0.149 ±0.093 


1.98 


0.722 ± 0.563 


0.154 ±0.289 


0.865 ± 0.297 


0.446 ±0.178 


tagged 


1.16 


-0.373 ± 0.293 


-0.363 ± 0.302 


-0.372 ±0.191 


0.119 ±0.174 


pairs 


1.10 


-0.079 ±0.196 


0.282 ±0.122 



Table 4: Results for (^-){pt) for the three categories of events, both targets and hadron charges. 
Only statistical uncertainties are shown. 
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Figure 15: The value of (^-){pt) determined in the anti-tagged category for protons (top) and 
deuterons (bottom) and positive (full points) and negative (open points) hadrons as a function of 
Pt- Also shown are the values for the tagged (squares) and pairs (triangle) category at their average 
respective pr- The uncertainties shown are statistical only. 



6.3 x dependence of ^ 

As there is no assumption- free method to determine the average x, various functional forms 



for ^-{x) with free parameters were investigated for extracting -^-(x) from Eqs. 4.5 and 4.6 



Assuming a functional form it is possible to convert the px dependence of the asymmetry 
into a value of -^-(%) at an average x. In contrast to the method described above, this 
method works for stronger x dependences of —2. 

For a given functional form and parameter set for ^-{x) Afj^px) can be calculated 



using Eq. 4.5. The best parameter set is obtained by minimizing the quantity 



X 



2 - (AA) T C A AA, 



(6.2) 



where AA is a vector containing the difference between the measured and the calculated 
asymmetries for each bin in pt 



AA = A r 



I f>BG aBG , nSIG A SIG\ 
\ a MC n A MC ) 



(6.3) 



The matrix C A in Eq. 6.2 is the covariance matrix including the statistical uncertainties 
of the data and MC asymmetries. 

A scan over an appropriately large parameter space is performed in order to find the 
parameters of the function describing ^-{x) that minimize \ 2 - Their covariance matrix 
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C F is determined from the distribution of probabilities at each scan point k: 



CS = — -^77 (6-4) 



k 



In this expression, 9^ is the value of parameter i at point k, while Q™ ax is the value of 
parameter i with the maximum probability. The probabilities can be evaluated from 
the x 2 cumulative distribution function. The advantage of this scan procedure is that it 
ensures finding the global minimum and enables the determination of the average x of the 
measurement using the extracted shape of 

This determination of — (x) can be done only for the anti-tagged category because of 
the necessity of having several bins in pr- In order to satisfy the fundamental requirement 
for — (x) to vanish at x = the functions are required to behave asymptotically as — (x) — > 

A//. 



x as x — > 0. In addition, lim — (x) — > 1 was required [78]. Omitting this constraint does 

ec-»l Q 



not significantly change the results. The small number of PT(beam) bins available limits the 
choice of the functional forms to those with no more than two free parameters. Several 
functional forms were studied, and the following two selected: 

fct. 1: x(l+pi(l - x) 2 ), 

fct. 2: x(l+pi(l-x) 2 +p 2 (l-^) 3 )- 
The parameters are restricted such that the LO positivity constraint: |-^(x)| < 1 is sat- 



isfied. Figure 16 compares the pt dependence of the measured asymmetry with the 
asymmetry calculated using the functional forms fitted to the measured asymmetries us- 
ing Eqs. and |4.5| . The \ 2 P er degree of freedom is large for both functions because 
of the discrepancy between the measured and calculated asymmetries in the highest px 
bin. No functional form was found that also accommodates the fourth data point within 
the statistical uncertainty. Systematic uncertainties of the Monte Carlo simulation (see 
Sect. |6.4| ) have not been used in this minimization; including them would reduce the x 2 
value significantly. 



Figure 17 shows the two functional forms of ~^-{x) an d their statistical uncertainties. 



The parameter value and uncertainties for fct. 1 are given in table |. The light shaded 
area represents the full x range spanned by the data, 0.07 < x < 0.7 (see Fig. |l3|). The 
dark shaded area represents the range of x spanned by preponderance of the data as seen 
in Fig. |l^. Although there are considerable differences in the ^-(x) functional forms over 
the full x range, the resulting Monte Carlo asymmetries are not very different, as can be 
seen from Fig. 16. From the behavior of the measured asymmetries together with the 
variation of the x-distribution (see Fig. [l3|) with PT(beam) it can be seen that any smooth 
function that describes the data leads to ^-{x) for x < 0.2 either small and positive or 
slightly negative, and significantly positive at larger x. However, no function with so few 
parameters is able to describe the sudden change of ^(x) at x 0.2 required to match 
the measured asymmetry in the largest PT(beam) bin. The average is determined using 
the resulting ^-(x) 
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Figure 16: Measured asymmetries with statistical uncertainties in four px bins for the anti-tagged 
category and a deuterium target, compared to calculated asymmetries using the two functions. 
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(6.5) 



where the sum is over all MC hadrons k in the px range 1 GeV < px < 2.5 GeV. This 
average determines the average (x) of the distribution probed by this measurement using 
the mean value theorem for integration, i.e., (x) is the value of x at which {-^) = ~^(( x ))- 

Figure 18 shows the total uncertainty (light gray band) of ^f-{x) vs - x evaluated with 
fct. 1 in the px range 1.0 GeV < px < 2.5 GeV and -^-((x)). The difference between 
fct. 1 and fct. 2 is assigned as an additional systematic uncertainty on the results from 



fct. 1 included in sys — models (see Sect. |6.4| ). The values of (x) determined from the two 
functions differ by only 0.007. 

The value for the gluon polarization extracted for the anti-tagged category from the 
deuterium target at (x) = 0.22 and a scale (fi 2 ) = 1.35 GeV 2 is 

^S((x), (/r 2 )) = 0.049 ± 0.034(stat) ± 0.010(sys-exp)l° ^(sys-models). 
The scale (fj 2 ) was determined by averaging over the scale of all signal MC events. The 
details on the systematic uncertainties are listed in table B. 
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Figure 17: Functional forms used with the values and statistical uncertainty bands from the fits. 
Light shaded area: the total x range spanned by the data (see Fig. |l3|); dark shaded area: the range 
in x where the preponderance of the data lies. 



6.4 Systematic uncertainties of 

6.4.1 Pythia, Jetset, and helicity-dependent (averaged) PDFs 

At present there is no Monte Carlo code available beyond leading order that models all 
subprocesses relevant at the kinematics of this experiment. Therefore this analysis is 



limited to leading order. As explained in Sect. 5.1, the Pythia model was significantly im- 
proved to better describe the Hermes helicity- averaged data over a wide kinematic range. 
The model contributions to the systematic uncertainty ( l sys — models') are determined 
by varying the parameters controlling the helicity- averaged and helicity-dependent PDFs, 
the Pythia subprocess cross sections and Jetset fragmentation process, and the low-pj- 
asymmetry. An individual uncertainty contribution is determined as the difference between 
with the standard setting and ^ obtained with the alternate setting. Related types of 
uncertainties are grouped in classes: 'parton distribution functions', 'Pythia parameters', 
'low-pr asymmetry', and 'fit function fct. 2'. All the individual and combined uncertainties 
are shown in table ||. 

For most types of uncertainties within a class, e.g., helicity-dependent nucleon PDFs, 
the uncertainty is conservatively estimated to be the maximum deviation appearing among 
the alternative models tested. Within a class these maximum differences are added in 
quadrature to form the 'Total' uncertainty for each class. The 'Total sys-models' uncer- 
tainty is obtained by adding those of all classes linearly, because of the complexity of 
correlations between them. 

Each of the classes investigated is motivated and discussed below. 
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Figure 18: The light gray band shows the total uncertainty of —f{ x ) vs - x with the statistical 
and total systematic uncertainty (see Table ||) added in quadrature. Note that the total systematic 
uncertainty contains a component accounting for the difference between fct. 1 and fct. 2. The point 
shown represents — {(x)) at (x) — 0.22. The inner error bar represents the statistical uncertainty 
and the outer the total uncertainty obtained by adding statistical and total systematic uncertainty 
in quadrature. 



• Parton distribution functions 
Spin- dependent nucleon PDFs 

The alternative parameterizations for the quark helicity distributions to GRSV, which 
were used are: GS-B [79], BB-06 [80], which includes the most recent g\ data from 
Refs. | 0, and the GRSV standard scenario [fTq] . The GS-B and BB-06 parame- 
terizations result in deviations, which are the second largest systematic uncertainty. 
They are of opposite sign and similar magnitude. The third alternative has a negli- 
gible effect. 

Spin- dependent photon PDFs 

Alternative parameterizations chosen are the maximum and minimum scenarios of 
GRS ]75|, [76| . The resulting deviations are of opposite sign and similar magnitude 
and make a significant contribution to the overall uncertainty. 
Spin-averaged nucleon PDFs 

Using the alternative parameterization GRV98 [74] for the spin-averaged quark and 
gluon distributions results in a small deviation. 
Spin-averaged photon PDFs 

The alternative parameterization is GRS [77], which results in a small deviation. 
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• Pythia parameters 

According to the discussion in Sect. 5.2 the following Pythia and Jetset parameters 
have been varied around their central values: PARP(91) and PARP(99), respectively 
the initial kj< of the partons in the nucleon and photon, are varied together; PARJ(21), 
which regulates p^ ag , where the upper and lower values correspond to an increase 
of xttd by 1 un it with respect to the standard setting |j6|; PARP(34), which is the 
multiplicative factor for the factorization and renormalization scales (PARP(34) p?). 
The scale in the calculation of the asymmetries was varied accordingly. In order 
to vary the subprocess fractions directly PARP(90) was varied. It is a parameter 

/ nPARP(90) 

regulating the cutoff pr min = PARP(81) ( parp(89) ) between direct and 

anomalous processes as well as soft and hard GVMD processes. 

The combined uncertainty of this class is comparable to that from the PDFs class, 
with PARJ(21) being the largest single contribution. 

• Low-pr asymmetry 

The logical alternative to A[ ow _ PT = D ■ A\ (which fits the HERMES low-p^ data) 
is to assume that at low pt any spin dependence is washed out, e.g., Ai ow _ PT = 0. 
Any such reduction of the asymmetry would only affect the lowest two px bins and 
could only increase 

• Fit function fct. 2 

For the x dependence, there is an additional class corresponding to functions with 
the shape of fct. 2. 



6.4.2 Experimental systematics 

The experimental systematic uncertainty is dominated by the fractional uncertainties in 
beam and target polarization, as shown in table [l]. They are added in quadrature and 
amount to 3.9% for the asymmetry and 20% for ^-((x)) from the deuterium target (shown 
in tables [| and ||) . 

Due to the rapid reversal of the target spin orientation (~ 90 s) the asymmetry ex- 
traction is independent of detector efficiency fluctuations. Possible false asymmetries due 
to the luminosity normalization are found to be negligible. 



6.5 Comparison to world data and models 

Only a few results obtained in leptoproduction exist on at present [26, 30, |3l], |32"|| . They 
were obtained from experiments with widely different kinematics and they have different 
scales fi 2 . Therefore, they cannot be easily compared. Nevertheless, for comparison the 
measurements are shown together at their respective (x) value, neglecting the Q 2 depen- 
dence of The experimental results shown in Fig. 19 are all obtained in leptoproduction, 
in LO analyses, although for different final states. The Hermes result is plotted with a 



horizontal bar indicating the half width at half maximum of the x distribution from Fig. 13 



Fit function fct. 1 is shown for the full x range spanned by the Hermes data (see Figs. 19 
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Table 5: Average kinematics and results for —(a;) and the parameter pi for fct. 1 with their 
statistical and systematic uncertainties, from deuteron data for anti-tagged events and for the full 
range in pr- For most types of uncertainties within a class, e.g., helicity-dependent nucleon PDFs, 
the uncertainty is conservatively estimated to be the maximum deviation appearing among the 
alternative models tested. Within a class (separated box), these maximum differences are added in 
quadrature to form the 'Total' uncertainty for each class, and referred to as 'Total PDFs' and 'Total 
PYTHIA'. These components, the 'low-p-r asymmetry' and the 'fit function fct. 2' uncertainties are 
added linearly to form the 'Total sys-models' uncertainty. 



and [l7|). The statistical precision of the Hermes result is the best currently available. 
The published Compass result for high-pf hadron pairs in the region Q 2 < 1 GeV 2 [31] 
has almost twice the statistical uncertainty. Concerning the systematic uncertainty the 
Hermes result is solidly based on varying many parameters of the well-tuned Pythia sim- 
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ulation and a comparison of results from several event categories and targets. The other 
results on ^ [26, 3C, [3^] are characterized by much larger statistical uncertainties. The 
earlier HERMES result of p^] is omitted, because the model used in this paper neglects 



important underlying subprocesses contributing to the signal and background asymmetries 
in the kinematic region used to extract ^f-{x) and no systematic uncertainty for the model 
used was determined. Altogether, the presently available experimental information from 
leptoproduction clearly indicates small values of over the covered x range. This con- 
clusion is consistent with the most recent results from polarized pp collisions from Phenix 
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Figure 19: The gluon polarization ^-((x)) from Hermes extracted with fct. 1 ((a;) = 0.22, (y 2 ) = 
1.35 GeV 2 ) compared to the ones from Compass [|l], M_M (low Q 2 : y 2 = 3 GeV 2 , high Q 2 : 
H 2 = 2.4 GeV 2 , open charm: y 2 = 13 GeV 2 ) and SMC ||0i 2 = 3.6 GeV 2 ) including statistical 
uncertainties (inner error bars) and total uncertainties (outer error bars). The x region of the data 
is indicated by the horizontal bars. Fit function fct. 1 is shown over the full x range spanned by 
the Hermes data. Also shown are a sample of curves from NLO pQCD fits DSSV, and BB-09) at 



1.5 GeV 2 . For clarity only the central values are shown. 



Also shown in Fig. 19 is ^{ x i calculated from two NLO pQCD fits, obtained as the 
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ratio of the helicity-dependent PDFs (DSSV U, and BB-09 [|8l| ) to the helicity-averaged 
PDFs (MRST02 @ and DR-09 ||], respectively). The BB-09 NLO pQCD fit is based 
on the inclusive DIS world data set. The DSSV NLO pQCD fit includes the world data 
on inclusive, semi-inclusive DIS and polarised proton proton scattering. The existing data 
obtained in leptoproduction on -^-{x) disfavor large magnitudes of the gluon polarization 
over the measured x range, in agreement with NLO-QCD fits by DSSV ||. 



7. Summary 

The gluon polarization in the nucleon has been determined by measuring the longitudi- 
nal double-spin asymmetry of high-py electro-produced single inclusive hadrons at Her- 
mes on a deuterium target. The value of —2- has been extracted using the measured 
asymmetries along with the subprocesses fractions, asymmetries and kinematics for the 
signal and background processes calculated using the leading-order Pythia Monte Carlo 
code. The value of is obtained from the product of gluon polarization and subpro- 
cess asymmetries summed over a wide range in x. The values of (^)(px) an d ^( x ) 
were determined independently. The systematic uncertainty was evaluated by varying 
Pythia parameters and models of background asymmetries. The final result for the av- 
erage gluon polarization in the pt range 1.0 GeV < p T < 2.5 GeV is ^f{{x), (fi 2 )) = 
0M9±0m4(stat)±0M0(sys-exp)±^ll(sys-models) at (x) = 0.22 and (/i 2 ) = 1.35 GeV 2 . 
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A. Cross section and asymmetries 



The leading-order formulas for helicity-dependent and helicity-averaged cross sections for 
longitudinally polarized virtual photons and partons for PGF, QCDC, DIS [84] and some 
QCD 2 — s> 2 subprocess are shown below. These are integrated over the azimuthal angle 
between the positron scattering plane and the production plane. The hard subprocess 
asymmetry (a) is given by 



a(s,i,n 2 ) = Aa/(2a), 
Act = <7 + — <r~ , 
a = (a+ + a~)/2, 



(A.l) 
(A.2) 
(A.3) 



with (+) denoting that both partons have the same helicity and (— ) the opposite helicity. 
The charge of the struck quark is given by e„ in units of the elementary charge. The 
common factor C is 4?r a e m 



DIS : 

d 2 a eq ^ q 
dudQ 2 

d 2 a eq ^ q 
dudQ 2 
d 2 Aa ec '^ q 

dudQ 2 



TCe g x, 



TC2e 2 q x. 



(A.4) 
(A.5) 
(A.6) 



For PGF and QCDC, there are helicity-averaged transverse and longitudinal, as well 
as a helicity-dependent transverse hard cross sections, note in the following d&i = ^| 



PGF : 



d&l 9 ^ qq 



dal 9 ^ qq 



dAa^ 9 ^ qq 



C- 



4vr (Q 2 + s) 2 (s + Q 2 ) 2 ' 



a s e 2 
C q 



1 



4vr (Q 2 + s) 2 
1 



2vr [Q 2 + sf 



Q 4 + s 2 u 2 + t 
(s + Q 2 ) 2 ut 

Q 2 



n 



s u z + r 



s + Q 2 



id 



(A.7) 
(A.8) 
(A.9) 



where s = (p q +Pq) 2 , t = (q — p q ) 2 , u = {q — pq) 2 and q,p q ,Pq are the 4-momenta of photon, 
final quark and anti-quark, as shown in Fig. |l|b. 



QCDC : 



da\ q ^ qg 



dal q ^ qg 



2a s e 2 
C q 



AQ 2 u 



3^ {Q 2 + s) 2 (s + Q 



2\2 • 



2a s e 2 
C q 



1 



3vr (Q 2 + s) 
1 



2uQ 2 



Q A + u 2 



(s + Q 2 ) 2 



Aa s e 2 
C q 



3tt {Q 2 + s) 2 



2(Q 2 -u) Q 4 + u 
Q 2 + s si 



st 

C.2 



(A.10) 
(A.ll) 
(A.12) 
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Reaction 



qg 
qg 
gg 
gg 



qg 
qg 
qq 
gg 



q a qb ->• q a qb 
q a qb ->• g'cgd 

qq -> fff 



da/di 



(* 2 + u 2 )[i 

f,2 i i2 



9 /o tfi 
2^ ~ F 



-3-0 

9sui 

-3-1 

3 P+u 2 
8 s 2 



f2 



St 



2~ 



3M ' 



* 2 -2 

32 P+ti 2 _ 8£ 
27 fij 3_ 



7T- 



dAa/di 



2(u 2 



2(u 2 -S 2 )[ 



9f/.s 



3 P+u 2 u 2 +P 



*1 
*1 



3 tit 



ut * 



9(-3 + 2§ 7 + F 



8 \s 



s -2 72 



9 L t 2 
8 

dadoed 



Sabi 

2 



-s 2 7 2s 2 yi 



[— ^ac^MTJ- — 8 r .dSab^-?r — h 



+ SabSbd 



W 

22 , 



3st 



27 fit 



+ 



t 2J 
16 P+u 2 
3 ~ P 



Table 6: Subprocess differential cross sections (A)a ab ~* cd for parton-parton interactions. The 

2 

common factor of has been omitted. 



For hard QCD 2 — > 2 processes the formulas [85] relevant in this analysis are given in 
Tab. §. 

For lepton scattering, the helicity-dependent and helicity- averaged cross sections are 
given by 



- DT T - T 



dvdQ 2 dsdi 

d A a eq^e'ff 

dvdQ 2 dsdt 



dsdt 



r T ( 



d 2 a, 



T 



2~7*9->-// 



dsdt 



+ e 



dsdt 



(A.13) 
(A.14) 



where / stands for a quark, an antiquark or a gluon in the final state and D is the virtual- 
photon depolarization factor 



D(y,Q 2 



y[(l + 7 V2) (2 - y) - 2y 2 m 2 JQ 2 



y 2 (1 - 2m 2 /Q 2 ) (1 + ^ + 2 (l + R)(l- y - 7 2 y 2/ 4) > 



(A.15) 



7 = Q jv , R = (Jl/pt for the subprocess, and is the transverse photon flux factor 



2vrQ 2 zy 



y 2 (1 - 2m 2 /Q 2 ^ 



2(l-y- 7 V/4) 
1+7 2 



and 



1 + - (1 - 2m 2 /Q 2 ) - 



y 2 + i 2 y 2 

_y _ 7 2 y 2/ 4 



(A.16) 



(A.17) 
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B. Tuned Pythia parameters 



MSEL=2 


MSTP(13)=2 


MSTP(17)=6 


MSTP(20)= 


=4 


MSTP(38)=4 


MSTP(61)=0 


MSTP(71)=0 


MSTP(81)= 


=0 


MSTP(92)=4 


MSTP(101)=1 


MSTP(121)=1 






PARP(2)=7 


PARP(18)=0.17 


PARP(91)=0.40 


PARP(93)= 


=2 


PARP(99)=0.40 


PARP(102)=0.5 


PARP(103)=0.5 


PARP(104) 


=0.3 


PARP(111)=0 


PARP(121)=2 


PARP(161)=3.00 


PARP(162) 


=24.6 


PARP(163)=18.8 


PARP(165)=0.477 


PARP(166)=0.676 






PARJ(1)=0.029 


PARJ(2)=0.283 


PARJ(3)=1.20 


PARJ(21)= 


0.40 


PARJ(41)=1.94 


PARJ(42)=0.544 


PARJ(45)=1.05 






MSTJ(12)=1 


MSTJ(45)=4 








MSTU(112)=4 


MSTU(113)=4 


MSTU(114)=4 






CKIN(1)=1.0 


CKIN(65)=1. • 10~ 9 


CKIN(66)=100. 







Table 7: The PYTHIA parameters, tuned to HERMES data, which are different from the default 
settings that can be found in Ref. [jf9| . 
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C. Systematic Uncertainties 





PT bin (GeV) 


1.0-1.2 


1.2-1.5 


1.5-1.8 


1.8-2.5 


1.0-2.5 


Kinematics 


(pt) (GeV) 


1.11 


1.30 


1.60 


1.90 


1.20 


Values 


PT 


&g/g 


0.017 


0.033 


0.026 


0.619 


0.055 


dependence 


5{Ag/g)(stat) 


0.067 


0.046 


0.073 


0.154 


0.033 


Systematic Uncertainties 


Category 


Model 




spin-dependent 


BB-06 


0.038 


0.019 


0.026 


0.089 


0.027 


Nucleon PDF 


GS-B 


-0.030 


-0.031 


-0.045 


-0.104 


-0.037 




GRSV-val 


0.007 


0.005 


0.008 


0.001 


0.006 


spin-dependent 


GRV(max) 


0.027 


0.019 


0.014 


0.060 


0.025 


Photon PDF 


GRV(min) 


-0.019 


-0.016 


-0.013 


-0.069 


-0.020 


Nucleon PDF 


GRV98 


0.020 


-0.008 


-0.031 


-0.144 


-0.006 


Photon PDF 


GRS 


0.029 


0.001 


-0.018 


0.100 


0.004 


Pythia 


PARP(90)=0.14 


-0.006 


-0.017 


0.003 


-0.076 


-0.016 


Parameters 


PARP(90)=0.18 


0.023 


-0.006 


-0.018 


0.014 


0.008 




PARP(91/99)=0.36 


0.002 


-0.002 


0.002 


-0.058 


-0.001 




PARP(91/99)=0.44 


0.004 


-0.009 


0.002 


-0.011 


-0.002 




PARJ(21)=0.38 


0.036 


0.018 


0.031 


0.110 


0.021 




PARJ(21)=0.42 


-0.023 


-0.040 


-0.032 


-0.072 


-0.034 




PARP(34)=0.5 


0.017 


-0.006 


-0.024 


-0.187 


-0.014 




PARP(34)=2.0 


-0.013 


0.012 


0.012 


0.210 


0.012 


1ow-pt asymmetry 


0.108 


0.037 


0.009 


0.005 


0.046 


Total PDFs ± 


0.058 


0.038 


0.059 


0.215 


0.045 


Total Pythia± 


0.047 


0.038 


0.047 


0.258 


0.033 


Total sys-Models + 


0.212 


0.113 


0.116 


0.473 


0.124 


Total sys-Models — 


0.105 


0.076 


0.107 


0.477 


0.078 


Experimental Systematic 


0.003 


0.006 


0.009 


0.002 


0.011 



Table 8: Average kinematics and results for (—){pt) with their statistical and systematic uncer- 
tainties, from deuteron data for anti-tagged events shown for the four bins and the full range in px- 
For most types of uncertainties within a class, e.g., spin-dependent nucleon PDFs, the uncertainty 
is conservatively estimated to be the maximum deviation appearing among the alternative models 
tested. Within a class (separated box), these maximum differences are added in quadrature to 
form the 'Total' uncertainty for each class, and referred to as 'Total PDFs' and 'Total PYTHIA'. 
These components and the i A low ~ PT ' uncertainties are added linearly into the 'Total sys-Models' 
uncertainty. 
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D. Hadron Asymmetries 



PT- 


bin 


(pt) 










deuterium 
















hydrogen 










(GeV) 


(GeV) 






(h + ) 




(h-) 






(h + ) 






(h-) 


n 


00- 


-0 1 5 


0.12 





nnn9 


4 


n 
u 


001 1 


n nn98 

U.UUiO 


4 


n 
u 


0011 





01 ^9 


4 


n 
u 


001 8 

\J\J 1U 


n 


ni dq 


4 


n 

u 


0020 


n 




u.ou 


0.23 


-0 


DDI 3 


4 


n 
u 


0004 


n 009^ 

U.UUiO 


4 


n 
u 


0004 





D1 D7 


4 


n 
u 


0007 

UuU I 


n 


UUOl 


4 


n 

u 


0008 


n 


^0- 
ou 


-fl 4^ 


0.37 


-0 


orn9 


4- 


n 


0004 

\J\J\J'~± 


n nn49 


4 


n 

u 


0004 





DDQ8 


4 


n 

u 


0007 


n 


uuou 


4 


n 

u 


0008 


n 




n fin 


0.52 


-0 


nn9Q 


4- 


n 

u 


0006 




-1- 

ZL 


n 
u 


0006 





m m 


4- 


n 
u 


0009 


n 


uuyj 


4- 

ZIZ 


n 
u 


0011 


n 


60- 


-0.75 


0.67 


-0 


0038 


4 


o 


0008 


-0.0037 


4 


o 


0009 





0127 


4 


o 


0014 


n 


0119 


± 


o 


0017 


o 


75- 


0.90 


0.81 





0005 


± 





0012 


-0.0027 


± 





0015 





0146 


± 





0021 


o 


0105 


± 





0026 





90- 


-1.05 


0.96 


-0 


0003 


± 





0019 


-0.0007 


± 





0023 





0166 


± 





0033 





0088 


± 





0043 


1 


05- 


-1.20 


1.11 





0069 


± 





0032 


-0.0033 


± 





0038 





0351 


± 





0055 





0091 


± 





0071 


1 


20- 


-1.35 


1.26 





0150 


± 





0054 


-0.0021 


± 





0063 





0563 


± 





0094 





0167 


± 





0118 


1 


35- 


-1.50 


1.41 





0174 


± 





0091 


0.0062 


± 





0104 





0487 


± 





0157 





0035 


± 





0197 


1 


50- 


-1.65 


1.56 





0429 


± 





0148 


-0.0017 


± 





0172 





0886 


± 





0256 


-0 


0759 


± 





0327 


1 


65- 


-1.80 


1.71 





0719 


± 





0238 


-0.0001 


± 





0277 





1317 


± 





0412 


-0 


0398 


± 





0530 


1 


80- 


-2.00 


1.88 


-0 


0075 


± 





0342 


-0.0027 


± 





0410 





1605 


± 





0596 





0428 


± 





0776 


2 


00- 


-2.50 


2.16 





0377 


± 





0463 


-0.0908 


± 





0572 





1071 


± 





0807 





0575 


± 





1128 


2 


50- 


-5.00 


3.04 


-0 


0071 


± 





0817 


0.1201 


± 





0960 





0112 


± 





1405 


-0 


0102 


± 





1546 



Table 9: Anti-tagged inclusive hadrons: measured longitudinal double-spin asymmetry for positive 
and negative hadrons on a deuterium (hydrogen) target. The uncertainty shown is statistical 
only. There is an additional overall normalization uncertainty of 3.9% (5.2%), all other systematic 
uncertainties are negligible. 
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PT- 


bin 


(pt) 






deut 


srium 










hydrogen 






(GeV) 


(GeV) 


A meas (h+) 


A meas (h-) 


A meas (h+) 


A meas {h-) 


0.00- 


-0.15 


0.10 


0.0229 


± 


0.0045 


0.0159 


± 


0.0050 


0.0896 


± 


0.0073 


0.0643 


± 


0.0089 


0.15- 


-0.30 


0.23 


0.0201 


± 


0.0029 


0.0120 


± 


0.0033 


0.0850 


± 


0.0048 


0.0605 


± 


0.0058 


0.30- 


-0.45 


0.37 


0.0188 


± 


0.0031 


0.0134 


± 


0.0036 


0.0807 


± 


0.0053 


0.0479 


± 


0.0064 


0.45- 


-0.60 


0.52 


0.0186 


± 


0.0040 


0.0134 


± 


0.0048 


0.0789 


± 


0.0068 


0.0639 


± 


0.0085 


0.60- 


-0.75 


0.67 


0.0211 


± 


0.0057 


0.0066 


± 


0.0068 


0.0661 


± 


0.0095 


0.0494 


± 


0.0122 


0.75- 


-1.00 


0.85 


0.0104 


± 


0.0072 


0.0106 


± 


0.0088 


0.0722 


± 


0.0122 


0.0683 


± 


0.0159 


1.00- 


-1.30 


1.11 


0.0465 


± 


0.0143 


0.0367 


± 


0.0176 


0.1055 


± 


0.0237 


0.1318 


± 


0.0314 


1.30- 


-1.60 


1.40 


0.0660 


± 


0.0364 


-0.0586 


± 


0.0431 


0.1410 


± 


0.0607 


0.0481 


± 


0.0813 


1.60- 


-2.00 


1.72 


0.0165 


± 


0.0903 


-0.0390 


± 


0.1086 


0.0627 


± 


0.1501 


-0.0229 


± 


0.2018 


2.00- 


-3.50 


2.18 


0.1534 


± 


0.3059 


0.3929 


± 


0.3798 


-1.5868 


± 


0.9548 


1.4778 


± 


1.4182 



Table 10: Tagged inclusive hadrons: measured longitudinal double-spin asymmetry for positive 
and negative hadrons on a deuterium (hydrogen) target. The uncertainty shown is statistical 
only. There is an additional overall normalization uncertainty of 3.9% (5.2%), all other systematic 
uncertainties are negligible. 
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A meas (hh) 




(GeV 2 ) 


hyc 


rogen 


deuterium 


1.00 


0.023 ± 


0.004 


-0.005 


± 0.003 


1.20 


0.018 ± 


0.006 


-0.005 


± 0.003 


1.40 


0.016 ± 


0.008 


-0.000 


± 0.004 


1.60 


0.031 ± 


0.010 


0.001 


± 0.006 


1.80 


0.047 ± 


0.013 


-0.006 


± 0.007 


2.00 


0.041 ± 


0.016 


-0.010 


± 0.009 


2.50 


0.084 ± 


0.028 


-0.028 


± 0.015 


3.00 


-0.001 ± 


0.045 


-0.041 


± 0.025 


4.00 


0.042 ± 


0.108 


-0.080 


± 0.057 



Table 11: Inclusive hadron pairs: measured longitudinal double-spin asymmetry for proton and 
deuterium targets. The uncertainty shown is statistical only. There is an additional overall nor- 
malization uncertainty of 5.2% (3.9%) for hydrogen (deuterium), all other systematic uncertainties 
are negligible. 
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